Kagome ice state in the dipolar spin ice Dy 2 Ti 2 07 
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We have investigated the kagome ice behavior of the dipolar spin-ice compound Dy2Ti207 in 
magnetic field along a [1 f 1] direction using neutron scattering and Monte Carlo simulations. The 
spin correlations show that the kagome ice behavior predicted for the nearest-neighbor (NN) in- 
teracting model, where the field induces dimensional reduction and spins are frustrated in each 
two-dimensional kagome lattice, occurs in the dipole interacting system. The spins freeze at low 
temperatures within the macroscopically degenerate ground states of the NN model. 

PACS numbers: 75.10.Hk, 05.50.+q, 75.25,+z 

Geometrically frustrated spin systems have been inves- 
tigated for decades because of their fascinating magnetic 
properties 0, 0, 0, 0, 0] ■ An Ising model on a cubic py- 
rochlore lattice consisting of spins parallel to local (111) 
easy axes was shown to exhibit macroscopic degeneracy 
in the ground states hi the same way as the pro- 
ton disorder in water ice does 0]. The local easy axis, 
ferromagnetic nearest-neighbor (NN) exchange interac- 
tion, and the geometry of corner sharing tetrahedra give 
rise to frustration in this system. Since the local "two- 
in and two-out" spin configuration of the ground states 
simulates the "ice rule" in water ice, the Ising model is 
referred to as "spin ice" 0] . 

Frustrated behavior attributable to the NN ferromag- 
netic spin-ice model was discovered in pyrochlore ox- 
ides, e.g. R2T12O7, R2S112O7 (R=Dy, Ho) 00,003, 
where the observed residual entropy agrees with the Paul- 
ing estimate for water ice 0- However, spins in these 
oxides interact dominantly by the dipolar interaction. 
Hence, they are represented by another Ising model with 
a dominant dipolar interaction, referred to as dipolar 
spin ice It was a puzzle why the dipolar spin ice 

shows the spin ice behavior predicted for the NN model. 
This was solved by Monte Carlo (MC) simulation studies 
[ill IT2I , showing that the dipolar interaction brings 
about an effective NN ferromagnetic coupling. In ad- 
dition, it was shown that the long-range dipolar inter- 
action only slightly lifts the degeneracy of the ground 
state manifold of the NN spin-ice model, and spins freeze 
within this manifold at low temperatures and never reach 
the true ground state [l^ . Il3| . The simulations also re- 
produce experimental observations of the dipolar spin-ice 
compounds More recently, an elegant analyt- 

ical understanding has been exploited using an argument 
of projective equivalence (lif . 

By applying magnetic fields to the NN spin-ice model 
the macroscopic degeneracy of the ground states can be 



partly or fully lifted depending on the field direction and 
mag nitude LZIUJ- For a field parallel to a [111] direction, 
along which the pyrochlore lattice is alternating stacking 
of kagome and triangular layers [see Fig. [fla)] , the field 
with intermediate strengths partly lifts the degeneracy 
and simultaneously induces decoupling of spins between 
the layers. In these ground states, spins on the trian- 
gular layers are parallel to the field, and spins on the 
kagome layers retain macroscopic degeneracy preserving 
the "two-in and two-out" configuration at each tetra- 
hedron [see Fig. mb)] 0. The ground states on each 
kagome layer can be mapped onto those of an exactly 
solvable dimer model in two dimensions 

Emms At 

low temperatures, by neglecting excited states with an 
energy gap, one can expect ice-like behavior of the spins 
on the kagome lattice, referred to as "kagome ice" |l7jj . 

For the dipolar spin ice DyaTioCv. experiments with 
field along [111] were performed [ll El El IH • It was 
shown that the magnetization plateau is consistent with 
the predictions of the kagome ice behavior for the NN 
model [rH E^]. However, the residual entropy plateau, 
another important point for proving the kagome ice, has 
not been definitely shown by the specific heat measure- 
ments 0, 0, [25| because of experimental difficulties 
[2l| . In addition, effects of the long range nature of the 
dipolar interaction have not been clarified for the kagome 
ice state. In this work, we performed neutron scattering 
experiments on Dy2Ti20y with field along [111] in the 
conjectured kagome ice regime, aiming at microscopic 
demonstration of the kagome ice state. We compared the 
neutron data with MC simulations based on the dipo- 
lar spin- ice model [ill ll^j ] . The simulations reproduce 
the main characteristics of the kagome ice: the residual 
entropy, the spin correlations close to the wave vector 
(§,— 1,0), and the field-induced dimensional reduction, 
being consistent with observations. 

Most of neutron scattering experiments were per- 
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formed on the triple-axis spectrometer GPTAS installed 
at JRR3, Japan Atomic Energy Agency. A single crystal 
sample of Dy 2 Ti 2 07 used in this work was 22 x 3.1 x 0.58 
mm 3 in size. Its long direction is parallel to a [111] di- 
rection, which ensures negligibly small demagnetization 
effect. The sample was mounted in a dilution refrigera- 
tor so as to measure the scattering plane perpendicular 
to the [111] direction. A scan along the [111] direction 
was carried out on the lifting-arm two-axis diffractometer 
D23 CEA-CRG at Institute Laue-Langevin. All the data 
shown are corrected for background and absorption. The 
MC simulations were carried out using a supercomputer 
at Institute for Solid State Physics, University of Tokyo. 
A standard Metropolis algorithm with single-spin-flip dy- 
namics was used in our simulations. The long-range dipo- 
lar interaction was handle d by a method |2 ll| equivalent 
to the Ewald method [ill H3f . We carried out the sim- 
ulations up to 38400 spins with up to 6 x 10 5 MC steps 
per spin. 

We adopted the Hamiltonian used in Refs. [ill lifjj 
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where Sf represents the spin vector parallel to the lo- 
cal (111) easy axis with \Sf \ = 1 at the sublattice site 
a in the unit cell of the fee lattice site i. The first term 
is the Zeeman interaction between the spins possessing 
the effective moment fi c g ~ 10 /ie and the magnetic field 
H. The second and third terms are the exchange and 
dipolar interactions, respectively. For the NN ferromag- 
netic spin-ice model there is only one exchange coupling 
constant J\ > and D = 0. 

It is well established that the low temperature behav- 
ior of Dy2Ti2C>7 can be approximately reproduced by a 
single-spin -flip MC simulation based on the dipolar spin- 
ice model [HI [l^, [Ri [25| except a few observations. For 
this model the dipolar interaction with D = 1.41 K and 
the Zeeman interaction are the dominant terms. The 
weaker exchange interactions consist of the NN antifer- 
romagnetic interaction with J\ = —3.72 K, and smaller 
second and third neighbor couplings with J2 = 0.1 and 
J3 = —0.03 K, respectively. The very weak J3 becomes 
important for determining the ground state within the 
degenerate states when the field is along [112 ] , because 
larger interactions balance each other |25j, |26( . Another 
balancing occurs for the [111] field close to the spin flop 
transition [see Fig. Hp)], which allowed us to determine 
J2 by a fit of the T dependence of the specific heat 
C(H = 0.75T) [see Fig.QJd)] to the simulation 0. 

The kagome ice behavior in Dv2Ti2Q7 was suggested 
by the residual entropy measurement 0, for fields 
close to H = 0.5 T as shown in Fig.^c). We calculated 
the residual entropy using MC simulation of the dipolar 
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FIG. 1: (color online) (a) The pyrochlore lattice is alternat- 
ing stacking of kagome and triangular layers along a [111] 
direction, (b) A spin configuration of the kagome ice state is 
shown, (c) The observed |2ljl and calculated residual entropies 
are plotted as a function of magnetic field. The calculated en- 
tropy for the NN model is S ( T = 0.1K), which is temper- 
ature dependent in the peaked region. Dotted lines represent 
the precise values of the spin ice and kagome ice states for the 
NN model [Tsl ITsll . (d) Temperature dependence of observed 
PH and calculated specific heat at H = 0, 0.5, 0.75 T. 



spin-ice model, where C(H)/T is integrated over 0.2 < 
T < 200 K. The MC calculations are shown in Fig. [He) 
and (d) together with those usin g th e NN spin-ice model 
with ferromagnetic J\ — 4.5 K 20]. The computation 
agrees with the experiment around H ~ 0.5 T, and shows 
more clearly the plateau behavior than the experiment 
|2l) . The low temperature entropy S(T = 0.1K) of the 
NN spin-ice model [joj . showing a peaked structure at 
a higher field (H ~ 1.35 T), indicates that the discon- 
tinuous decrease of the residual entropy at H ~ 0.9 T, 
accompanying the spin flop transition [23| . is brought 
about by the long-range dipolar interaction. 

The plateau of the residual entropy around H ~ 0.5 
T of the dipolar spin-ice model is slightly smaller than 
that of the NN model. This suggests that the kagome-ice 
ground-state manifold of the NN model is weakly lifted by 
the dipolar interaction. For the origin of the kagome ice 
behavior of Dy2Ti20y one naturally expects that spins 
freeze within these nearly degenerate ground states ow- 
ing to the high energy barriers for breaking the "two-in 
and two-out" ice rule Il2l. This mechanism can be mi- 
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FIG. 2: (color) (I) HT phase diagram of Dy2Ti207 obtained from the magnetization and the specific heat measurements |2lll23|. 
The kagome ice state is observed in the blue region. The red solid line represents the first-order spin-flop transition. (II, III, 
IV) Neutron intensity patterns in the scattering plane perpendicular to the [Iff] direction are shown for several temperatures 
and fields denoted by crosses in (I). In (II-a)~(II-f) experimental results for Dy2Ti207 are shown. Calculations using the MC 
simulation based on the dipolar and NN spin ice models are shown in (III-a)~(III-d) and (IV-a)~(IV-d), respectively. 



croscopically investigated by observing spin correlations 
using neutron scattering techniques. 

In Fig. Olll) we show neutron intensity patterns in 
the scattering plane perpendicular to the [111] direction 
at several points in the HT plane denoted by crosses 
in Fig. Ell). The low temperature data at T — 0.3 K 
show definite differences between the spin ice state at 
H = [(II-a)], the kagome ice state at H = 0.5,0.8 
T [(II-b),(II-e)], and the fully aligned state at H = 1 
T [(II-f)], where very weak scattering intensity implies 
that very small spin fluctuations are left above the spin 
flop transition. We note that the intensity peak observed 
for the kagome ice state [(II-b),(II-e)] at the wave vector 
Qo = (§, — §7 0), denoted by black circles, is a character- 
istic of the kagome ice state. This spin correlation was 
found in the exact ground states of the NN spin- ice model 


In Fig. 0111) an d (IV) we show scattering intensities 
calculated by MC simulations based on the dipolar and 
NN spin-ice models, respectively. The simulations for the 
spin ice state [(III- a),( IV-a)1, which are consistent with 
previous results |14ll27| . confirm that the agreement with 
experiments [(II-a)] is improved by taking account of the 



dipolar interaction. In contrast to this, the simulations 
for the kagome ice state [(III-b),(IV-b)] show almost the 
same pattern to the naked eye, and agree well with the ex- 
periment [(II-b)]. As noticed previously, the field induces 
decoupling of the spins between the kagome layers for the 
NN model. This field-induced dimensional-reduction is 
confirmed by no Q dependence of the calculated inten- 
sity pattern [(IV-b)] along the [111] direction. For the 
dipolar model, we also observed little Q dependence of 
the calculated intensity pattern [(III-b)] along the [111] 
direction, and hence conclude that the spin correlations 
show quasi-two dimensional character in the kagome ice 
state. As temperature is increased at the typical field 
H = 0.5 T of the kagome ice state, the experimental in- 
tensity patterns [(II-c),(II-d)] are consistent with the sim- 
ulations [(III-c,d),(IV-c,d)]. From these results we con- 
clude that the degeneracy of the kagome ice ground states 
of the NN spin-ice model is almost identically preserved 
for the dipolar spin-ice model and for Dy2Ti2C>7, provid- 
ing microscopic evidence for the kagome ice behavior of 
this compound. 

In order to study the spin correlations around Qo a t 
H = 0.5 T, we performed scans along the (h, —h, 0) line 
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work (ill Hif explaining the equivalence of the NN and 
dipolar models for the spin ice state, there may be a sim- 
ilar enlightening argument also for the kagome ice state. 

In summary, we have investigated the kagome ice be- 
havior and field-induced dimensional reduction of the 
dipolar spin ice Dy2Ti2C>7 in magnetic field along the 
[111] direction using neutron scattering. The observed 
spin correlations were analyzed by MC simulations based 
on the dipolar and NN spin- ice models. For fields around 
H ~ 0.5 T spins freeze within the degenerate kagome-ice 
ground states of the NN model, which are weakly lifted by 
the dipolar interaction, and consequently exhibit kagome 
ice behavior. 



FIG. 3: (color online) Q-scans along the (h, —h, 0) line 
through Qo at H = 0.5 T: (a) experimental and (b) MC sim- 
ulation results based on the dipolar model. Curves are fits to 
Lorentzian functions, (c) Temperature dependence of the in- 
verse correlation length £ _1 obtained by the fits. The dotted 
lines are guides to eyes. The solid line represents the expo- 
nential T dependence exp(— A/2fc_gT) described in the text, 
(d) Q-scan along the [111] direction from (211) [Fig.^II-b)]. 



and plot the results in Fig. [3J where the MC calculations 
of the dipolar model are also shown. One can see from 
Fig.[3Ja) that the spin correlations of Dy2Ti2C>7 cease to 
develop around T = 0.5 K. This disagrees with the theory 
of the NN model , where it was argued that the cor- 
relation length diverges exponentially as exp(A/2/csT) 
(T — > 0), where A is a gap energy, in the thermal equi- 
librium state. The non divergence of the spin correlations 
for Dy2Ti207 can be ascribed to spin freezing due to the 
barrier for breaking the ice rule. As discussed in Ref. 01 
for the spin ice state, the single-spin-flip dynamics of the 
MC simulation become non-ergodic at low temperatures, 
and hence, the simulated system does not reach a ther- 
mal equilibrium state within practical MC steps, as well 
as Dy2Ti2 07 does not within experimental time scales. 
In fact, the simulated system also show freezing around 
T = 0.7 K, which does not necessarily agree with the ex- 
periment because of difference of the spin-flip dynamics 
between the simulated system and Dy2Ti20y. To observe 
the two-dimensional spin-correlations in the kagome ice 
state, we performed a Q-scan along the [111] direction 
from (211) [see Fig.HJII-b)]. The result together with the 
MC simulations of the two models are shown in Fig.jSfd). 
These experimental data, being consistent with the sim- 
ulations within the experimental error, suggest the two 
dimensional nature of the kagome ice state in Dy2Ti20y. 

Finally we note that a surprising fact revealed by the 
present numerical study is the high similarity of the 
spin correlations, shown in Fig.|2lII-b)^(III-d) and (IV- 
b)~(IV-d), between the NN and dipolar models espe- 
cially for the kagome ice state. In view of the theoretical 
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